Collective Excitations in Quantum Hall Liquid Crystals: Single-Mode Approximation Calculations 
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A variety of recent experiments probing the low-temperature transport properties of quantum Hall systems 
have suggested an interpretation in terms of liquid crystalline mesophases dubbed quantum Hall liquid crystals. 
The single mode approximation (SMA) has been a useful tool for the determination of the excitation spectra of 
various systems such as phonons in ^He and in the fractional quantum Hall effect. In this paper we calculate (via 
the SMA) the spectrum of collective excitations in a quantum Hall liquid crystal by considering nematic, tetratic, 
and hexatic generalizations of Laughlin's trial wave function having two-, four- and six-fold broken rotational 
symmetry, respectively. In the limit of zero wavevector q the dispersion of these modes is singular, with a gap 
that is dependent on the direction along which q = is approached for nematic and tetratic liquid crystalline 
states, but remains regular in the hexatic state, as permitted by the fourth order wavevector dependence of the 
(projected) oscillator strength and static structure factor. 

PACS numbers: 73.43.-f, 73.20.Mf, 64.70.Md. 52.27.Aj 



I. INTRODUCTION 

For more than two decades two-dimensional electron sys- 
tems (2DES), and in particular the quantum Hall effect 
(QHE)ii2i^ have been constant sources of complex and unex- 
pected behavior, perhaps with no equal in the realms of con- 
densed matter physics. The unique combination of extremely 
high mobilities (/i ~ 10^ mA^s) in GaAs/Al^^Gai-a^As het- 
erostructures, low temperatures (T < 100 mk), enhancement 
of interactions due to the reduced dimensionality, and relative 
quenching of the kinetic energy in strong magnetic fields due 
to Landau level (LL) quantization, has allowed the emergence 
of complex and striking behavior due to subtle correlation ef- 
fects. In fact, from the discovery of the integer and fractional 
QHE's in the early eighties^ (leading to the award of two No- 
bel prizes), to the existence of novel fractionally chargecP- and 
composite particles^ and many newer interesting many-body 
phenomena, the QHE has been a consistently active and ex- 
citing area of research. 

The physics of 2DES in partially filled LL's is inherently 
complex due to the high degeneracy of the "unperturbed" 
ground state (i.e. without the Coulomb interaction). The first 
successful theoretical approach to this system was proposed 
shortly after the discovery of the FQHE^ by Laughlin, who 
proposed his famous trial- wavefunction- 
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to describe states at filling factor v = 1/m, with m an odd 
integer. Here Zj = i 



iyj is the position of the j*^ electron posed in the pasti 



Later, Jain proposed^ the elegant conceptual framework 
of the composite fermions (CF) which unifies all hierarchies 
of integer and fractional QHE's, along with the intermediate 
states in between QH plateaus: an even number 2p of vor- 
tices is attached to each electron, also lowering the Coulomb 
energy, and (in a "mean field" approximation) the "magnetic 
fluxes" associated with the vortices lead to a reduction of the 
effective magnetic field: B* = B — 2p^one (^o = h/e is 
the flux quantum and Ue is the electron density), resulting in 
an effective filling factor v* = — 2pv). It is easy to see 
that for the strongest FQHE states z/* is an integer (e.g. for 
V = l/Z and p = Iv* = 1), leading to the interpretation of 
the FQHE of electrons as a simple IQHE of the CF's; whereas 
the intermediate regions (e.g. u = 1/2, 1/4) correspond to 
z^* ^ oo, i.e. a vanishing effective magnetic field: B* = 0. 

Since 1999, magnetotransport experiments have uncovered 
a variety of surprising results at low temperatures (T ^ 100 
mK), for example: extreme anisotropics^ and apparent com- 
petition between different ordered phases^ in the intermediate 
regions between quantum Hall plateaus at high LL's, the melt- 
ing transition between a Wigner crystal (WC) at u 1/7^ 
and the microwave conductivity evidence of structural phase 
transitions in partially filled LL's.^ A large body of evidence 
corresponding to seemingly distinct phenomena may be partly 
understood in terms of a simple conceptual point: that the 
many-body states involved have an intrinsic crystalline or 
liquid crystalline orderpiiSiiiiiSii^iiii^ii^ be it smectic^i^ 
nematic tetratic or hexaticw^ 

Generalizations of Laughlin wavefunction [Eq. ([U] with 
discrete broken rotational symmetry (BRS) have been pro- 
candidates for nematic or hex- 



in the complex plane and we work in units of the magnetic 
length /q = [Ti/eBY^'^. The "goodness" of Laughlin state 
originates in that the nodal hyper- surfaces on which the many- 
body wavefunction vanishes coincide with the ones in which 
the particles are in contact. Whereas the vanishing of the 
wavefunction when two (spin aligned) electrons are in con- 
tact is required by Fermi statistics, Laughlin state has multiple 
nodes at those points, thus reducing the Coulomb repulsion. 



atic statesiSi^ in order to understand anisotropic transport ob- 
served in the intermediate regions,^ or the melting of the WC 
at z/ = 1/7.^ In fact, the motivation for these states arises from 
the fact that it is generally expected that melting in 2D may 
occur through a topological Koszterlitz-Thouless-type (KT) 
transition. For 2D melting, the reliable Kosterlitz-Thouless- 
Halperin-Nelson- Young (KTHNY) theory ^^-^^ predicts that, 
in fact, an intermediate liquid crystalline phase will exist be- 
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tween a solid and a liquid phase, which will exhibit no trans- 
lational order, and only a quasi-long-range order for the orien- 
tational order below the KT transition temperature. These ar- 
guments were used by Wexler and Dorsey^^ to calculate qual- 
itatively correct anisotropic-isotropic transition temperatures 
for the quantum Hall liquid crystal in the transitional regions 
at high LL's.^ 

In this paper we consider the spectrum of collective excita- 
tions for a family of liquid crystalline states in a partially filled 
LL. These states are generated so as to satisfy the follow- 
ing criteria which we consider reasonable for understanding 
the dynamics: (i) states must obey Fermi statistics, i.e., they 
must be odd under the exchange of any pair of electrons; (ii) 
the states must be translationally invariant (far enough from 
boundaries); (iii) there must be a broken rotational symme- 
try belonging to the proper symmetry group (i.e., C2 for a 
nematic, C4 for a tetratic, and Cq for a hexatic; additional 
symmetries are possible in principle, e.g. with a Cio symme- 
try, we have not explored such possibilities); (iv) states and 
excitations must reside entirely in the LLL to avoid the large 
cyclotron energy cost hcOc. Note, as we will show later, that 
the properties of any excited LL may be readily obtained from 
the properties of the LLL. 

States that satisfy the abovementioned requirements have 
been proposed and studied in detaiUSJiiS for filling factors 
u = 1/3, 1/5, and 1/7. These are found by splitting the "ex- 
tra" vortices of the Laughlin (or other CF) states around the 
electron, while obeying the required symmetries: 
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where the complex numbers are distributed in pairs of op- 
posite value in the complex plane (to satisfy Fermi statistics). 
For the states with the highest discrete symmetry at each fill- 
ing factor we may take 



an = a e 



i 27r(/i-l)/2p 



/iG{l,2,...,2p}, 



(3) 



and without loss of generality a can be taken to be real. 
The wavef unction in Eq. (|2|l represents a homogeneous liquid 
crystalline state at = l/(2p + 1), is anti- symmetric, lies en- 
tirely in the LLL, and is smoothly connected to the isotropic 
Laughlin state for a = 0. Figure [l] depicts the nodal distri- 
bution for various states of Eq. A remarkable feature of 
these states is that they posses an underlying charge density 
wave (CDW), but these CDWs are melted by fluctuations, and 
overall the system is translationally invarianti^ 



II. THE SINGLE MODE APPROXIMATION 

To calculate the excitation spectrum we use the single mode 
approximation (SMA),^^'^^ which reliably provides the first 
moment (mean) of the energy of the excitations (for a given 
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FIG. 1: Nodal distribution for Zij = Zi — zj for a quantum Hall 

nematic at z/ = 1/3, tetratic at zy = 1/5, and hexatic at zy = 1/7. 

wavevector k) that are coupled to the ground state by means of 
the density operator.^^'^^'^^ The SMA was first used by Feyn- 
man in 1953 to accurately calculate the spectrum of phonons 
in superfluid ^He.^^ The essence of the method originates on 
the assumption that the ground state of a system of bosons 
has a scarcity of long- wavelength excitations. Under those 
circumstances, the variational wavef unction for an excitation 
corresponding to a density-wave can be written as 



/)k(ri, ...,yn)=N -^/^PkV^o(ri, . . . ,rAr) , 



(4) 



where pk 



2^7 = 1 ' 



ikr, 



is the density operator, and ?/^o is 



the many-body ground state (which is, in fact, unknown for 
^He). Note that this trial state automatically builds in the fa- 
vorable correlations of the ground state. The energy of this 
excited state, A (k) = ((/)k|i^ — £^o|</>k)/(</>k|^k), can be sim- 
ply evaluated: 



A(k) 



5(k) 



(5) 



In the last term the numerator is the "oscillator strength" and 
takes on the universal value /(k) = Ti^k'^ /2m, and 5'(k) 
is the static structure factor, which is directly measurable by 
means of neutron scattering (it is here, in fact, where the He- 
He correlations of the ground state are "encoded"). Using 
the experimental results for 6'(k), Feynman could calculate 
a spectrum of remarkable quality, showing the phonon-like 
spectrum at small wavevector s, and a roton minimum at wave- 
vectors comparable to the inverse interatomic distance 



III. THE SINGLE MODE APPROXIMATION IN THE 
QUANTUM HALL EFFECT 

The applicability of the SMA to fermion systems is also 
well established for two- and three-dimensional systems in the 
absence of magnetic fields, giving a good approximation for 
the plasmons at long wavelength. For 2DES in presence of 
a magnetic field it correctly gives the zero-wavevector mag- 
netoplasmon at cjc = eB/rUe, a result that is guaranteed by 
Kohn's theorempi which states that the dipolar excitation is 
saturated by the cyclotron mode (this results in the modes of 
interest — the intra-LL excitations — having quadrupolar char- 
acter, i.e. with an oscillator strength / (x q^). 

For excitations fully contained within a single LL, the cy- 
clotron mode is not of primary interest. In 1985 Girvin, Mac- 
Donald an Platzman (GMP) proposed an ingenious ansatz for 
projected excited states 
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FIG. 2: (Color online) Pair correlation function ^(r) for a = 1/3 nematic, a = 1/5 tetratic, and a = 1/7 hexatic. In all cases 
5'(r) = (9[r^/^] for small r and ^(oo) ^ 1. The fact that ^(r) ~ 1 over a large area is a guarantee that the calculations are able to reproduce 
the thermodynamic limit. 



iV^q) = Pq IV^O) , 

where pq is the projected density operatoriS^ 

= ^{0,m\e-''^-''\0,m)al^,ao,n 



(6) 
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AT N 



where |0, m) correspond to single-particle states in the lowest 
LL and angular momentum m, and aj ^ is the particle creator 
operator for such state. As in Feynman's ansatz [Eq. ©I, Eq. 
^ preserves the favorable correlations of the ground state. 
The exclusion of inter-LL excitations eliminates the problem 
with the saturation of the dipolar mode. The excited states 
have a compelling description, in first quantized formi2^ 
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(8) 



X tlj{zi, . . . , Zj-i,Zj - iq, Zj^i, . . . , zn) , 

which corresponds to shifting each electron by x q and 
superimposing these N configurations with an amplitude 

^-iq*Zj/2^ 

Similarly to Eq. (|5ll, the excitation spectrum can be readily 
obtained 



(2iV)-i(V>o|[pt,[g,Pq]]|^o) _ /(q) 

" 5(q) ■ 



(9) 



The projected oscillator strength comes from the non- 
commutation of the projected density operator with terms in 
the potential energy part of the Hamiltonian also projected 
onto the LLL: 



Since [pk, Pq] = (e^*^/^ - e^^*/^)Pk+q, we findSS 



/(q) 



2e-l^l'/2^sin2(^)5(k)x 
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(11) 



For its part, the projected static structure factor S'(q) can be 
calculated fromi22i2^ 

5(q) = 5(q) - (1 - e-«'/2) ^ (^2) 
where S'(q) is the unprojected static structure factor: 

5(q)-l = po / cPre-^^-'[g{v)-l], (13) 



(14) 



the Fourier transform of the pair correlation function 
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which is obtainable from the ground state via, e.g., Monte 
Carlo (MC) simulations JSiliSi In our case, we considered 
BRS states [Eq. Q] corresponding to a = 1/3 nematic, a 
v = 1/5 tetratic, and slu = 1/7 hexatic. Figure |2| depicts the 
pair correlation function for various states. In all cases, corre- 
lation functions and SMA excitation spectra were computed 
for numerous a's. The angle dependence significantly in- 
creases the burden in the MC simulations since the full angle- 
dependent ^(r) is needed, rather than the considerably simpler 
angle-averaged_^(r) of the isotropic systems. The accurate 
calculation of /(q), with its angular-dependent exponentially 
large factors required high quality 5'(q) and hence ^(r). To 
put things in perspective, to achieve ca. 1% accuracy in ^(r), 
(9 [10, 000] counts were accumulated for each 0.01 x 0.01 (in 
units of lo) r-box in the original histograms (the small boxes 
were necessary to have precision in the Fourier transform at 
high wavevectors). This required, for each filling factor u and 
anisotropy generating a, to run approximately 4-8x10^ MC 
steps in systems of = 200-400 electrons taking 100-200 
cpuxdays of com puta tion for each (z/, a) in a 2 GHz Athlon 
cluster (see Ref. Iisll for further details). A relatively large 
A^e is required so that the simulations are able to reproduce a 
system in the thermodynamic limit. The fact that ^(r) :^ 1 
over a large area is a guarantee of that achievement (Fig.|2|). 

Figure [3] shows the projected static structure factors S{q) 
for a z^ = 1/3 nematic, v = 1/5 tetratic,s and sl v = 1/7 
hexatic. From S{q), the oscillator strength /(q) is computed 
using Eq. fTTl . Analysis of S{q) and / (q) shows that both are 
for small q.^^ This restriction on the the small wavevec- 
tor behavior originates in Kohn's theorem,^^ as all of the O [q'^] 
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V = 1/3, a= 1.5 



V = 1/5, a = 2.0 





FIG. 3: 
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(Color online) Projected static structure factors S (q) for a. v ■ 
■ 0[q^] for small q, and 5'(oo) ^ 0. 




pieces in the unprojected parts are saturated by the (uninter- 
esting) inter-LL excitations at fiuOc- 

Figure |4| presents some of our results for the excitation 
spectra in the lowest LL. The results are consistent with those 
obtained by GMP for the isotropic v = l/Z and 1/5 FQHE 
cases^^ which were qualitatively confirmed experimentally.^^ 
They show that the collective excitation spectrum remains 
gaped, albeit with a deep magnetoroton, for modes coupled 
to the ground state via the density operator. Not surprisingly, 
BRS states have significant anisotropy in their spectra. How- 
ever, it is noteworthy that for the nematic and tetratic cases the 
spectrum is singular, with an angle dependence on the exci- 
tation energy A(q) as q ^ 0. By contrast, the hexatic liquid 
crystal has a regular spectrum in the long wavelength limit. 
The apparent disparity have, of course, to do with the differ- 
ent rotational symmetries of the different states: as A = 
and both numerator and denominator are 0[q^] for small q, 
there is no possibility of generating a Cq symmetric form with 
terms that can only depend on q^, Qy, and q^qy. This can also 
be understood from the point of view of an effective elasticity 
theory for 2DES (valid in the long wavelength limit, g ^ 0) 
the elasticity tensor being of the fourth rank is not compatible 
with a six-fold rotational symmetry. This elastic interpretation 
will be published elsewhere. 

The presence of this singular spectrum suggests that mi- 
crowave conductivity experiments2i2£ may be able to discern 
such structures as a signature of, e.g., the quantum Hall ne- 
matic suggested by magnetotransport experiments.^'^ 

Another interesting feature (Fig.|3 is that the magnetoroton 
minima and the gap at the origin appears to collapse at some 
finite anisotropy factor a. The appearance of elastic modes 
that go soft may be a precursor of the appearance of charge 
density waves in the system. Spectra in higher LL's can be 
obtained by multiplying the projected density operator pq by 
^l(^^/2), where Ll{x) is the Laguerre polynomial and L 
corresponds to the desired LL^ The rightmost panel of Fig. 
|4] shows the modification of the SMA spectrum for the first 
excited LL. Numerical error due to the large LL{q^ /2) factors 
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= 1/3 nematic, a z/ = 1/5 tetratic, and a = 1/7 hexatic. In all cases, 

makes it difficult to get reliable results for higher LL's at this 
point. 

Theoretical predictions of which type of order is lowest in 
free energy (required to decide which state is favorable at fi- 
nite temperatures, as shown in Refs. |6 7]) are still incom- 
plete, however, as the entropy is likely to be dominated by 
gapless modes originating, e.g., from the Goldstone modes 
associated with the spontaneous breaking of the continuous 
rotational symmetry of the isotropic states. Generalization of 
the SMA to gapless rotational modes (see e.g. 15) will re- 
quire, however, 3 -body operators which demanding consider- 
ably higher computing capabilities^ 



IV. SUMMARY 

Summarizing: in connection with the recent variety of ex- 
perimental evidence supporting liquid crystalline phases in 
quantum Hall sy stems piSii^ we have calculated the collective 
excitation spectrum in the SMA approximation for a quantum 
Hall liquid crystal. We found that the spectrum of excitations 
coupled to the ground state by the density operator remains 
gapped, but develops a significant anisotropy, which in the 
case of the nematic and tetratic liquid crystals has a singular 
gap in the long wavelength limit. 
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